Atom-Molecule Coherence in a Bose-Einstein Condensate 



Elizabeth A. Donley, Neil R. Claussen, Sarah T. Thompson, and Carl E. Wieman 
JILA, University of Colorado and National Institute of Standards and Technology, Boulder, Colorado 80309-0440 

(Dated: February 1, 2008) 

Coherent coupling between atoms and molecules in a Bose-Einstein condensate (BEC) has been 
observed. Oscillations between atomic and molecular states were excited by sudden changes in 
the magnetic field near a Feshbach resonance and persisted for many periods of the oscillation. 
The oscillation frequency was measured over a large range of magnetic fields and is in excellent 
quantitative agreement with the energy difference between the colliding atom threshold energy and 
the energy of the bound molecular state. This agreement indicates that we have created a quantum 
superposition of atoms and diatomic molecules, which are chemically different species. 



There is considerable interest in extending the appli- 
cations of ultracold atoms to ultracold molecules. One 
route for producing a very cold and possibly Bosc- 
condensed sample of molecules is to create the molecules 
from an atomic BEC. Wynar et al. 1 created cold 87 Rb2 
molecules in a single ro-vibrational state of the ground- 
state potential from an 87 Rb BEC using a two-photon 
stimulated Raman transition. The authors could not 
probe the coherence properties of the molecules in that 
state, but the prospect of creating a superposition of 
atomic and molecular condensates initiated a flood of 
theoretical work on the subject 2-6 . Ultracold molecules 
have also recently been formed through photoassociation 
of a sodium BEC 7 . 

Utilizing the natural atom-molecule coupling that 
arises from a Feshbach resonance is an alternate route 
for producing ultracold molecules from an atomic BEC, 
and it is the route we have followed here. A Feshbach 
resonance is a scattering resonance for which the total 
energy of two colliding atoms is equal to the energy of 
a bound molecular state, and atom-molecule transitions 
can occur during a collision 8-12 . A schematic representa- 
tion of the potentials involved is shown in the inset of Fig. 
1A. For our 85 Rb resonance, BEC atoms in the F = 2, 
mp = —2 state collide on the open-channel threshold. 
F and mp are the total spin and spin-projection quan- 
tum numbers. The bound state in the closed channel 
differs in energy by an amount e from the open-channel 
threshold. The bound molecular wave function can be 
described as a sum of amplitudes of different hypcrfinc 
components {F, mp) having Mp = 771f,i + ~m-F,2 = — 4 13 . 
Because of their different spin configurations, the atoms 
and molecules generally have different magnetic moments 
and the difference depends on magnetic field. Thus e de- 
pends on magnetic field and the degree of atom-molecule 
coupling is magnetically tunable. The energy difference 
between the free atoms and the bound molecules is plot- 
ted in Fig. 1A. This behavior of the bound-state energy 
also causes a resonance in the scattering length, a, which 
is shown in Fig. IB. The scattering length characterizes 
the mean-field interaction energy of a BEC. 

When the magnetic field is tuned to values near the 
Feshbach resonance, theory predicts coherent oscillations 



between the atomic and molecular states, but there is 
significant disagreement on the conversion fraction and 
the coherence properties 16-21 . 

In experiments with a sodium BEC, Stenger et al. 22 ob- 
served that inelastic losses were dramatically enhanced 
when they ramped the magnetic field across the Fesh- 
bach resonance. We observed similar results for 85 Rb, 
but with lower rates 23 ' 24 . It is likely that the formation 
of molecules played a role in the loss, but there was no 
experimental evidence for the presence of molecules and 
the results followed a loss-rate dependence on time. More 
recently, we measured the time dependence of the loss 
in an 85 Rb BEC by applying controlled magnetic-field 
pulses toward but not across the Feshbach resonance 25 . 
We observed the surprising result that under some con- 
ditions, shorter, more rapid pulses actually led to more 
loss than longer, slower pulses that spent more time near 
the resonance. The time dependence of the loss was sug- 
gestive of a nonadiabatic mixing of states, with the only 
states within a reasonable energy range being the normal 
atomic BEC state and the nearby bound molecular state. 

In this work we show that much of the loss is likely due 
to the coherent mixing of atomic and molecular states. 
To create a superposition and probe its coherence, we 
applied two short magnetic-field pulses toward the Fes- 
hbach resonance, separated by a "free evolution" time 
during which the magnetic field was held at a constant 
value some distance from the resonance. We measured 
the number of atoms in the condensate as a function of 
time between the two pulses for various values of the 
steady-state magnetic field between the pulses. We ob- 
served dramatic oscillations in the number of atoms re- 
maining in the atomic BEC at frequencies corresponding 
to the energy splitting between the molecular and the 
atomic states. 

Experimental methods. The apparatus has been 
described in detail elsewhere 23,25 . We first created 85 Rb 
condensates 23 typically containing 16,500 atoms, with 
fewer than 1,000 uncondensed thermal atoms. The initial 
number Ni n u fluctuated from shot to shot by ~500 atoms 
(~ 3% number noise). After producing the condensate at 
a field of ~162 G, we ramped the magnetic field adiabat- 
ically to ^166 G, corresponding to an initial scattering 
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FIG. 1: Feshbach resonance bound-state energy and scatter- 
ing length. (A) Energy splitting versus magnetic field. The 
resonance is centered at ~155 G. The solid curve is a theo- 
retical estimate of the energy found with a coupled-channels 
calculation 14 , and the dotted line indicates e = 0. The in- 
set schematically shows the collision channels involved in the 
resonance, e depends on magnetic field because the atoms 
and molecules have different magnetic moments and thus the 
potentials have different Zeeman shifts. (B) Scattering length 
versus magnetic field for fields above the Feshbach resonance. 



length dinit — 10 ao, where ao = 0.053 nm. The spa- 
tial distribution of the atoms was Gaussian with a peak 
atom density of no = 5.4 x 10 13 cm~ 3 , and the trap fre- 
quencies were (17.4x17.4x6.8) Hz. After preparing the 
condensate we applied a selected fast magnetic-field pulse 
sequence by sending an appropriate time-dependent cur- 
rent through an auxiliary magnetic-field coil 25 . A typ- 
ical pulse sequence is shown in Fig. 2. It is composed 
of two nearly identical short trapezoidal pulses separated 
by a region of constant (but adjustable) magnetic field. 
Upon completion of the fast-pulse sequence in Fig. 2, 
we ramped the magnetic field from ~166 G to ~157 G 
in 5 ms and held at that field for an additional 7 ms 
to allow the repulsive mean-field energy to expand the 
condensate. Then we turned off the magnetic trap and 
used destructive absorption imaging 12.8 ms later to ob- 
serve the atomic condensate and measure the number 
of remaining atoms 26 . This detection scheme was nei- 
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FIG. 2: Magnetic field pulse shape. Fields shown for pulses 
#1 and #2 correspond to scattering lengths of ~2500 ao, and 
the free precession field B evo i ve corresponds to a scattering 
length of ~570 ao. The dashed line indicates the position 
of the Feshbach resonance. In the text, we refer to the free 
precession time as t evo i ve . The rise/fall time for all of the 
pulses that we used was 14 /is. 



ther sensitive to atoms with kinetic energies larger than 
~ 2 /iK nor to atoms in off-resonant molecular states. 
We determined the value of the magnetic field between 
the pulses, B evo i ve , by measuring the resonance frequency 
for transitions from the F = 2, mp = — 2 to the F = 2, 
rriF = —1 spin state by applying a 10 fis RF pulse to a 
trapped cloud of atoms 25 . 

Burst and remnant oscillations. As we observed 
for single pulses toward the Feshbach resonance 25 , there 
were two distinct components of atoms observed in the 
absorption images and a third "missing" component that 
we could not detect. One of the observed components was 
a cold remnant BEC which was not noticeably heated or 
excited by the fast-pulse sequence, while the other com- 
ponent was a relatively hot (~150 nK) "burst" of atoms 
that remained magnetically trapped during the BEC ex- 
pansion time. Using a variational approach 27 to model 
the mean-field expansion that we applied to the BEC 
remnant to measure its number, N remnant , we found that 
we should impart < 3 nK worth of energy to the remnant 
before imaging. This estimate agrees well with the ex- 
pansion velocity that we observed after the trap turn-off. 
Thus the remnant BEC was nearly 50 x colder than the 
burst. 

The missing component contained atoms that were in 
the initial sample but were not detected after the trap 
turn-off. To find the number of atoms in the remnant 
BEC and the number of burst atoms, we allowed the mag- 
netic trap to "focus" the burst cloud before imaging 28 . 
A typical image is shown in Fig. 3. We fit the focussed 
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FIG. 3: An absorption image taken after the fast magnetic- 
field pulse sequence and the mean-field expansion. The color- 
bar indicates the optical density. The horizontal and vertical 
directions coincide with the axial and radial axes of the trap, 
respectively. The dimensions of the image are 366 x 52 /im. 
The BEC remnant is the roughly spherical cloud at the cen- 
ter, while the burst atoms are focussed into a thin line along 
the axial direction. Note the dramatic difference between the 
two spatial distributions, owing to the large difference in their 
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burst (which had a much larger spatial extent than the 
remnant) with a two-dimensional Gaussian surface, ex- 
cluding the central region of the image that contained 
the remnant. This fit yielded the number of burst atoms, 
Nburst ■ Subtracting this fit from the image and perform- 
ing a pixel- by-pixel sum of the central region of the image 
gave N rernnan t . 

Nremnant versus t evo i ve is plotted in Fig. 4 for two dif- 
ferent values of B evo i ve . The number clearly oscillates. 
Changing the value of B evo i ve affected the oscillation fre- 
quency dramatically (note the change in scale from Fig. 
4A to Fig. 4B). After only pulse #1 and the subsequent 
constant field but with no pulse #2, N remnant showed no 
variation except for a slow decay consistent with the loss 
rate expected for a single pulse to that field 25 . 

We have taken data similar to those in Fig. 4 for a va- 
riety of different B evo i ve values. As in Fig. 4, we fit each 
curve to the function y = yo+A exp (—t/r) sin (2nvt + (p) 
to find the oscillation frequency v and decay time con- 
stant t. The measured frequencies are plotted versus 
Bevoive in Fig. 5 along with theoretical predictions for 
the bound-state energy relative to the atomic state. 

In the regime where the scattering length is much 
larger than the radius of the interatomic potential well, 
the bound state energy for an arbitrary attractive po- 
tential can be approximated by e = —h 2 /ma 2 29 . H is 
Plank's constant divided by 2n, m is the atomic mass, 
and a is the scattering length. The same equation re- 
lates the bound state energy to the effective scattering 
length, which is calculated from the Feshbach resonance 
parameters through the relation a = a^ g x (1 — B A B() ), 
where a^g is the background scattering length, A is the 
width of the Feshbach resonance, and Bq is the resonant 
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FIG. 4: N rem nant versus t evoive for n = 5.4 x 10 1 ' cm . 
(A) Be VO i V e = 159.69(4) G {devolve = 590 do). The data fit 
to a damped sine wave with a frequency of 207(2) kHz and a 
decay time of 46(21) us. The open squares near N rem nant = 
6000 indicate the number remaining versus time after only 
pulse #1 and tevoive at 159.69 G. (B) Bevoive = 157.60(4) G 
(devolve — 1390 do). Note the increase in time scale from (A). 
These data fit to an oscillation frequency of 23.9(12) kHz and 
a decay time of 82(38) /j,s. 



magnetic field 30 . 

The quantity \e\/h is plotted with no adjustable pa- 
rameters in Fig. 5. The measured oscillation frequen- 
cies are in excellent agreement with this simple model 
over the range of magnetic fields where the model is ex- 
pected to be valid. The theoretical results found with 
a much more sophisticated coupled-channels scattering 
calculation 14 in Fig. 5 are in excellent agreement with 
the data over the entire range. 

The fact that the oscillations occurred at exactly 
the frequency corresponding to the bound-state energy 
clearly indicates that we are creating a coherent superpo- 
sition of atoms and molecules with the sudden magnetic- 
field pulses 11,12 . Although we do not have a detailed 
understanding of how the field pulses couple atoms and 
molecules, by choosing the shapes of the perturbing 
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FIG. 5: Oscillation frequency versus magnetic field. The 
points are the measured frequencies. The solid line is 
the energy difference between the atom-atom threshold and 
the bound molecular state found by S. Kokkelmans with a 
coupled-channel scattering calculation. The dotted line is a 
plot of e/h. The inset is an expanded view of the lower- 
frequency data. The maximum frequency that we could mea- 
sure was limited only by timing jitter and finite resolution 
in the experiment. The magnetic-field measurements for the 
points with the smallest horizontal error bars were performed 
on the same days as the corresponding frequency measure- 
ments. The error bars for the points with larger field un- 
certainties were inflated by 100 mG to account for estimated 
day-to-day field drifts. 



pulses such that a single pulse results in roughly 50% 
loss, we observe high-contrast oscillations in the number 
of atoms in the atomic BEC. 

From the amplitude of the oscillations, one can put a 
lower bound on the number of molecules being created. 
Take, for example, the data in Fig. 4A. The amplitude 
of the atom oscillations was 1800(300) atoms. Assuming 
the fringes are coming from interference with molecules, 
there must be at least 1800/2 = 900(200) molecules on 
average. Assuming that the missing atoms are molecules 
that we fail to convert back into atoms gives an upper 
bound of 3200(100) molecules for the conditions of Fig. 
4A. 

The damping time for the oscillations, Tdecay, was more 
difficult to measure with high precision than the oscil- 
lation frequency. To within our measurement precision, 
Tdecay did not depend on B evo i ve , but our uncertainties in 
Tdecay were as large as 100% for some fields. We had the 
highest precision measurements for frequencies around 
200 kHz where the oscillation period was long compared 
to our experimental timing jitter but short compared 
to Tdecay At frequencies near 200 kHz, we measured 

Tdecay = 38(8) US for U = 1.3 X 10 14 Cm -3 and Tdecay = 

91(33) us when we decreased no to 1.1 x 10 13 cm -3 . 

Nburst also had interesting dependencies. For the con- 
ditions under which most of the data were collected, 



FIG. 6: Number versus t evo i ve for no = 1.1 x 10 13 cm"' ! . 
From bottom to top, the data are plots of Nburst (open 
circles), N remnan t (filled circles), and the total number of 
observed atoms (gray squares). Each data set was fit to 
a damped sine wave resulting in the displayed fits. Ninit 
= 17,100 is indicated by the flat dashed line. B evo i ve = 
159.84(2) G and the remnant data fit to an oscillation fre- 
quency of 196(1) kHz and Tdecay = 91(33) ^ts. To produce 
condensates with the lower density for these measurements, 
the fast-pulse sequence was applied from our evaporation field 
of 162.2 G and the amplitudes of pulses #1 and #2 were re- 
duced to ~7 G. 



the burst contained ^5000 atoms on average, which is 
~30% of Nina- N ourst depended on density and var- 
ied from one-half of the atoms lost from the condensate 
for our typical peak density of no = 5.4 x 10 13 cm" 3 
to nearly all of the atoms lost from the condensate for 
no = 1.1 X 10 13 cm" 3 . Nburst, N remnant , and total num- 
ber of atoms detected are plotted in Fig. 6 for B evo i ve 
= 159.84(2) G and n = 1.1 x 10 13 cm" 3 (5x lower den- 
sity than was used for the data shown in Fig. 4). All 
three components oscillated at the same frequency. The 
burst oscillation lagged behind the remnant oscillation by 
155(4)°. Since the relative phase shift is nearly 180°, the 
oscillation amplitude for the total number was smaller 
than cither the burst or the remnant oscillation ampli- 
tudes. The relative phase depended sensitively on the fall 
time of pulse #2. For example, when we increased the 
fall time from 11 fj,s to 159 /is, the burst oscillation then 
lagged behind the remnant oscillation by 68(7)° and the 
peak-to-peak amplitude of the total number oscillation 
was 5600(400). 

For the conditions of Fig. 6, Ninit — 17,100 exceeded 
the time-averaged total number of atoms counted after 
the pulse sequence by 8(3)% on average. For the higher- 
density measurements in Fig. 4, 39(4)% of the atoms 
were missing. Experiments with longer pulses #1 and ^2 
also had a higher fraction of missing atoms. For example, 
when we used 50 /xs pulses with no = 5.4 x 10 13 cm -3 , 
56(3)% were missing. 
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We have carried out double-pulse measurements with 
a variety of widths and amplitudes for pulses #1 and #2 
and a variety of different densities and initial magnetic 
fields. Although the oscillation frequency was unchanged, 
the phase, contrast, and damping of the oscillations did 
vary. The contrast was very sensitive to the pulse length, 
and was lower for longer pulses that created more missing 
atoms. Defining the contrast as the oscillation amplitude 
divided by the time averaged number of remnant atoms 
detected, we observed an optimum contrast of 0.42(7) 
for 15 fis pulses to 156.6(1) G. A single such pulse re- 
moved about half of the atoms from the BEC. When the 
pulse length was comparable to Tdecay, the contrast was 
reduced by about a factor of two. Under those condi- 
tions, ~3/4 of the atoms were lost after pulse #1. The 
phase was shifted, but we did not observe a change in the 
contrast when we varied the amplitudes of pulses #1 and 
#2 from B = 156.6(1) G (2,400 a ) to B = 155.1(1) G 
(24,000 ao). The contrast did depend on the intermediate 
level, however, and was reduced for B evo i ve values clos- 
est to the resonance, for which the magnetic-field jumps 
between B evo i ve and pulses #1 and #2 were shortest. 

We also looked for a temperature dependence of both 
the damping and the frequency at ^200 KHz and did not 
see any. The high-temperature data was much noisier 
than the data for pure condensates, due to unexplained 
enhanced noise in the number of thermal atoms after 
the magnetic-field pulse, but when the initial thermal 
fraction was increased from <5% to 30%, the data still fit 
to oscillations with frequency, amplitude, and damping 
consistent with what was observed with low temperature 
data. 

Conclusions and Outlook. Our interpretation of 
our observations is that the first magnetic-field pulse pro- 
vides a sufficiently rapid perturbation to result in nonadi- 
abatic mixing between atomic and molecular states. The 
superposition then evolves according to the energy dif- 
ference between the states, which is determined by the 
magnetic field during the free evolution stage, B evo i ve . 
The second pulse mixes atom and molecule states again, 
such that the final state of the system depends on the 
relative phase of atomic and molecular fields at the time 
of the second pulse. This is somewhat analogous to Ram- 
sey's method of separated oscillating fields 33 . Under very 
limited conditions (a near 1700 do), we could also observe 
Rabi-like oscillations with a single pulse towards the Fes- 
hbach resonance. This narrow window results from the 
conflicting needs for both strong coupling and condensate 
loss time 25 long compared to a Rabi oscillation period. 

After pulse #2, a fraction of the coherent molecular 
component is converted into the energetic but still spin- 
polarized burst atoms through a yet to be determined 
process. Another mystery is the missing atoms. Are they 
molecules that are not converted back into atoms and are 
not detected in the burst or the remnant signals? If so, 
why do we not see them as atoms after the field is turned 



off and the corresponding molecular state is no longer 
bound? Why are there fewer missing atoms for lower- 
density condensates and quicker pulses towards the Fesh- 
bach resonance? What is the actual conversion efficiency 
from atoms to molecules and how could we maximize it? 
Very near the Feshbach resonance, the molecular state 
has a magnetic moment nearly the same as that of the 
free atoms, and hence will remain magnetically trapped. 
A major remaining question concerns the nature of the 
molecules. Could they be considered a molecular BEC? 
Clearly there is much to be learned about this curious 
system. 
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